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ABSTRACT 

Linear state signal feedback is used to obtain exponential 
response from fourth order systems. Characteristic equation roots 
are selected to provide the desired exponential response with con- 
straint on initial conditions and system acceleration. A digital 
computer root locus program is developed to determine feedback 
coefficients in a manner which minimizes the possibility of oscilla- 
tory response in the presence of state sensor errors. The effect 
of noise on the state signal is investigated and a sample data fil- 
tering technique developed. A quasi-optimum time technique utilizing 
second order switching logic for initial control effort and linear 


feedback for terminal control is developed. 
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1.0 Introduction 

The purpose of this study was to investigate means of controlling 
fourth order systems to obtain various types of responses. The initial 
portion of the study deals with the use of limear state signal feedback 
to obtain exponential settling. Use of the uncoupled form is investi- 
gated for use in the feedback solution analysis. The effect of noise 
on the state signals for the linear feedback solution was investigated 
and a sample data filtering technique developed, A root locus program 
was developed to determine feedback coefficients in a manner which will 
minimize the possibility of oscillatory response in the presence of 
state sensor errors. 

The remainder of the study deals with a quasi-optimum time 
solution which incorporates switching logic. Full control effort is 
used initially and switching takes place at predetermined levels of 
state variable combinations, After switching, control is allowed to 
decay for a linear termination of the solution. The objective here 
is to obtain a near- optimum time response with no possibility of con- 
trol chatter. 

Section 2.0 of this report contains the linear feedback portion 
of the investigation. Section 3.0 contains the quasi-optimum investi- 
gation. All synthesis for the study was conducted on the school's 
Control Data Corporation 1604 digital computer using fortran. All 


programming will be referenced in the text and shown in the appendices. 


2.0 Linear Feedback for Exponential Settling 


The response of a system to a set of initial conditions may be con- 
trolled by feedback of the system state variables. That is, the forcing 
function of the system is made up of predetermined amounts of each of 
the system's state variables. Thus, the characteristic equation of the 
controlled system may be adjusted to obtain the desired response. 

If an exponential response of the system's position state variable 
is desired, the characteristic equation will have only negative real 
roots, with a dominant root that causes the desired exponential path 
after the decay of the non-dominant roots. The proper feedback co- 
efficients of each of the state variables may be calculated to give the 
desired characteristic equation. 

Example 

An aircraft landing flare is representative of a class of auto- 
matic control problems in which a system has initial conditions of each 
state variable and it is desired that the position state variable 
settle to zero in an exponential manner. In order to design a suit- 
able control system, a mathematical description of the aircraft longi- 
tudinal motion is required. Assuming constant air speed and a shallow 
glide angle leads to the short period equation of longitudinal motion 
2] These are written in terms of the following transfer function 
relating elevator position, d (radians), to aircraft pitch angle, 


p (radians); 
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where K = short period gain 
w = short period resonant frequency 
a = short period damping factor 
T = path time constant 
In order to complete the mathematical description of the aircraft, 
the transfer function relating altitude, h (feet), and pitch angle, 
p (radians), in terms of velocity V (feet per second) and path time 


constant T is 


h(s) = —————_ p(s) (2) 


s(Ts + 1) 


Combining (1) and (2) results in a transfer function relating 


altitude and elevator position: 





h(s) = id d(s 3 
- s“(s° + 2as + 7) - = 


Letting the system forcing function u be equal to KVd, 2 signal 


flow diagram of the system is shown in Figure l. 





Figure 1. Block Diagram of Aircraft Longitudinal Motion with 
State Variable Feedback 
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The characteristic equation for the system of Figure 1 is 
4 5 2\.2 
= ( 
s) + (b), + 2a)s~ + (b, +w)s + bs + b, = 0 (), ) 
It is now possible to solve for the feedback coefficients 


b,, b., b., ana b), such that the characteristic equation will have 


1? as 3 

the desired roots. Suppose, for a particular aircraft 

a = 0.5 and w = 1.0, and the desired closed loop characteristic equa-~ 
tion roots are s =-0.18, -1.0, -1.0, and -5.0; giving the charac- 


teristic equation 


s* + 7.1862 + 12.268" + 6.988 + 0.9 = 0. (5) 


Equating (4) and (5) will then give the desired feedback co- 
efficients 


b, = 0.9, b, = 6.98, b, = 11.26, b= 6.18. (6) 


This will give a time response solution for the flare, after allowing 
a short decay time (t, + 6 seconds) for the non-dominant roots, of 
approximately 

h(t) = h(e,) et (7) 
2.1 Use of the Uncoupled Form 

In order to more closely examine sources of instability in the 

above type of problem, and to determine the proper variable for a root 
locus study, it is useful to solve for the system's uncoupled form [3]. 


Rewriting (3) in the time domain gives 


(aye 2 eh(c) Ro h(t) mae) (8) 
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In order to rewrite (8) in matrix form let 


x,(t) = h(t) (9) 
rot) = b(t) 
x,(t) = h(t) 
h(t ) 


x(t) 
u(t) = Kvd(t) 


This leads to the following matrix form.for the system; 


0 ie Oo 0 O 
s-f8 8 hohe [he 0 
O 0 -wW -2a 1 
which is defined as 
x = Fx + Du (11) 


The system can now be transformed into the uncoupled form by 


defining a new variable y as 


y, = u(s)/s° (12) 
Y¥, = u(s)/s 

y, = u(s)/(s° + 2as + w°) 

y, = u(s)s/(s° + 2as + w°) 


After expanding x by partial fraction expansion it is noted that 


WEE -eb fas (hazaw=) Ju! Bays 
2 2 2 
nS L/w -2afwo -1fw |g fais) 
0 ) l 0 = 
0 On 0 1 
which is defined as 
x = Gy ° ; (14) 





Also, by solving for the inverse of G, the expression 


y= ox (15) 


ue 2A 1 O 
0 Ww ea 1 
[aoe Ow | ole (16) 
0 0 0 1 
In order to draw a signal flow diagram of the uncoupled system, 
(14) is substituted into (11) giving 
Gy = FGy + Du (17) 
; -1 
and multiplying by G 
; -l -1 
y=G FGy +G Du (18) 
Solving (18) gives the uncoupled form 
O 1 O O O 
: 0 O 0 OQ 1 | 
Le 0 0 0 eee Pe @ulae (19) 
O Oueen mee 1 


which is drawn as shown in Figure 2. Lines are added to denote the 


subsequent addition of state variable feedback loops. 











Figure 2. Block Diagram of Uncoupled Form of Aircraft 
Longitudinal Motion with State Variable Feedback 


The characteristic equation of the uncoupled system of Figure 2 is 


h 


Bua 2 
s° + (2atn,+n, )s + (wo4n) ¢2an,+n )s (20) 
2 o 
+ (2an, +w n,)s + nw = 0. 
Assuming the same aircraft as before with a = 0.5 and w = 1.0 gives 
3! + (l+n_+n )s> + (l+n,+n +n vee + (n,+n,)s +n, = 0 (21) 
ea 4 me 3 | Dae, 1 


Suppose, prior to beginning the flare, the aircraft is descend- 
ing a EIS glideslope with an airspeed of 170 knots and that the 
flare will begin at 100 feet. Nominal initial vertical velocity is 


approximately -20 feet per second. In order for the glide angle to 
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be approximately tangent to the flare path, a dominant characteristic 


equation root of -0.2 is selected. That is, the slowest phase plane 


eigenvector is placed through the nominal initial condition in the 


x1 Xp phase plane. Now the feedback coefficients for the control system 


can be calculated. Suppose characteristic equation roots are selected 


at s = -0.2, -1.0, -1.0, and -5.0. This gives the following desired 


characteristic equation: 


eu + 7.28? fe lene Fas fe Ih pa 78 


Equating (21) and (22) gives the following solution for the un- 


coupled state variable feedback coefficients 


n, = 4.2, aie OQ 


which can be written in matrix form as 


N = oe oan oy ee ae 8) ; 
The solution for the original state variable feedback coeffi- 
cients can be calculated by noting 


B = nev 


which yields 


Bee | 1.0 saa ee 6.2 | 


2.2 Use of a Root Locus Study 


In order to become familiar with the behavior of the roots of 
the characteristic equation (21), a fortran root locus program was 


written. Since the uncoupled characteristic equation contained n, 


as an element of each of the internal coefficients, n, was used as 


the variable for the root locus. As shown in Appendix I, it was 


(22) 


(23) 


(24 ) 
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found that (24) had all negative real roots for the dominant root vary- 
ing from -0.2 to -0.366. Thus the previous solution is on the boundary 
where a slight error in a feedback setting could cause two of the 

roots to leave the real axis and result in a oscillatory component 

in the system response. For this reason, it is better to solve for 

a dominant root of about -0.16 and then adjust n. to move the dom- 


a 


inant root back to -0.2. Doing this, as shown in Appendix I gives 


org Tee oe Orr || (27) 
[0.9 6.61 10.89 5.81 | (28) 


"a 
il 


and B 


These feedback coefficients give all real roots for a dominant 


root varying from -0.18 to -0.284 when n, is varied from 4.96 to 


2 
6.09. Figure 7 shows the system time response resulting from use of 
the state variable feedback coefficients of (28). 
2.3 Consideration of Acceleration Constraint 

The proper method for selecting the non-dominant root locations 
of the closed loop system characteristic equation has been ignored 
in the previous discussion. It will now be shown that system accelera- 
tion during exponential settling is a function of both initial con- 
ditions and the roots of the closed loop system characteristic eq- 
uation. Thus, if the most severe initial conditions that the system 
can be expected to be subjected to are predictable, an acceptable 
location for the non-dominant roots can be determined, 

As before, the dominant root should be selected to place the 
slow eigenvector through the nominal initial condition in the KyX5 


phase plane. As the non-dominant real roots of the closed loop 


system characteristic equation are moved to the left from the origin 
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in the s plane, the acceleration maximum during the early portion of 
the response is increased. Therefore, the non-dominant roots can be 
selected by considering the maximum amount of acceleration to which 
the system should be subjected. In the aircraft flare problem, the 
proper acceleration constraint would be determined by structural con- 
siderations as well as passenger comfort. 

In order to illustrate the effect of non-dominant root placement, 
a graphical solution of the acceleration on a fourth order system 
during exponential response was made. The closed loop fourth order 
system was generalized by lumping the feedback coefficients with the 
plant coefficients for each state variable. This resulted in a gen- 
eralized fourth order system as shown in Figure 3 with initial con- 


ditions on the state variables. 


hy (0) h,(0) h,(0) h, (0) 





—— 


Figure 3. Block Diagram of General Fourth Order System with 
Initial Conditions on the State Variables 


10 





Using signal flow techniques, the state variable transforms can 
be expressed in terms of the initial conditions on the state variable;. 
As is often the case for mechanical systems, the magnitude of maximum 
control effort will be determined by the acceleration state variable. 


The operational expression for the acceleration state is 


ah (0) +h (oie rc s“) - h,(0)(c_s+e,) - h,(O)se 
h(s) = 0 yes y 2 pag 1 1 (es) 


D 
S Cc, S Cc ._s + C.S§S + C 
jee aay 2 l 


3 


If a purely exponential response on the position state variable is 
desired, the roots of the closed loop system characteristic equation 


must be real. This can be expressed mathematically as 


sH4c,s3tc,8°4c,5+¢, =: (s+p, )(s+P, )(s+p,)(s+p, ) | (30) 


which leads to the following time expression for the system acceleration: 
2 
h(t) = €,04[h,(0)(P5P 3p), }#hy(0)(PpP4PP),+0 5Py, (31) 
h(0)(P#P#P, +h, (0)| ~ 
2 
e Pa) by (0)(P,PsP }483(0)(P,P5+P,P),+PP,, )# 
h,(0)(py+P #0), 4h, (0) + 
© 305] hy (0)( PoP, #83 (0)(P, Pat? 1+ Py) 
h(0)(p{ptp), +h, (0) | + 
e,,P; [hy (0)(PyPsP5 }4h(0)(p,Po+P,P,+PsP + 


b ,(0)(py+P #3 )+h, (0)| 


i 





where 


S Daatt 
al 


*1 = (pete s-P Mey) (32) 
Po“ Py P3 Py P), Py 


Curves may then be plotted for various initial conditions and 
characteristic equation roots. Figure 4 is a graphical solution 
obtained from the fortran program in Appendix II. It is a three 
dimensional plot of system acceleration resulting from various 
initial conditions on the position and velocity state as well as 
various non-dominant root locations. The dominant root location 
Psat -O.,2% 

2.4 Statistical Considerations 

This portion of the study was made to investigate the effect 
of noisy state variable signals when linear feedback is used. This 
would be the actual case in the previous aircraft flare example 
since a predictable amount of error due to the aircraft state sensors 
could be expected. Radio altimeter errors would be caused by thermal 
gaussian noise, terrain uneveness, and the effects of close proximity 
to the ground during the flare. The use of an altimeter as the 
prime sensor would necessitate differentiating three times in order 
to obtain four state variables. Because of the amount of noise that 
could be expected on the altimeter output, this would be a highly 
unlikely approach. Since an accelerometer would not be sensitive 
to terrain unevenness and proximity, a simulated accelerometer output 
was used as the prime system sensor for the computer synthesis. 


A simulated altimeter output was used as a secondary system sensor. 
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Figure 4 Acceleration of Fourth Order System with Dominant 
Characteristic Equation Root of 0.2 
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Sensor noise was simulated using a gaussian noise generator Le 
See Appendix III for the mathematical development and computer check- 
out of this noise generator. 

A digital computer was used to simulate the control problem. A 
fourth order Runge-Kutta numerical integration algorithm was used to 
simulate system motion. 

As in the previous example, it was desired that the system be con- 
trolled over an exponential flare. The plant of Figure 1 with 
a= 0.5 and w= 1.0 was used. The feedback coefficient of equation 
(28) were selected to give a dominant closed loop characteristic 
equation root at S = -0.2. 

In the computer synthesis, the simulated accelerometer output, 
with gaussian noise, is used to calculate the four system states, 
Smoothing is used to minimize the effects of the added noise. 
Periodically, the altitude state is updated with a simulated noisy 
altitude output which has also been smoothed. Updating is accomplish- 
ed by averaging the two altitudes. 

Smoothing of both the simulated acceleration and altitude is 
accomplished by calculating a least square error line over a vari- 
able number of past sensor outputs. See Appendix IV for a mathe- 
matical development of the smoothing methad. 

Figure 5 shows the time response of the system in a three 
dimensional graph. This multi-curve graph shows the effects of 
accelerometer noise variation on the response. Eleven Bivens are 
shown with accelerometer noise variance varying from zero to 0.01 


feet per second per second. Lines joining points of equal time 
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Fagune 5 Time Response of Fourth Order System with 
Noisy Sensors 
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(1 second, 2 seconds, etc.) are shown to give an illusion of depth to 
the graph. For each curve the altitude noise variance was 0.25 feet. 
Linear smoothing over three sample points, having a time duration of 
0.025 second apart, was used for both acceleration and altitude. 
These curves show that the terminal portion of the response becomes 
erratic when accelerometer variance is 0.003 feet per second per 
second or more. 

Figure 6 is similar to Figure 5 except that the third axis shows 
the effect of using a variable number of past sample points for linear 
smoothing. Eleven curves are shown with the number of sampling 
points varying from zero (no smoothing)to twelve, Acceleration 
variance was 0.005 feet per second per second and altitude variance 
was 0.25 feet. It is seen that as sampling points are increased to 
five or more, the data becomes sufficiently stale so that oscillation 
occurs. As the number of past sample points is increased even more, 


instability occurs. 
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3.0 Use of Second Order Switching Logic for Quasi- 

Optimal Response 

Minimum time solutions for specific fourth order systems can be 
found in the literature [ - These solutions are derived through 
the use of Pontryagin s maximum principle. 

At this time, a general solution for the fourth order minimum 
time problem has not been developed because of the high degree of 
complexity involved. This portion of the investigation was conducted 
to investigate the application of the well developed second order 
minimum time solutions to the fourth order problem i6] . 

The time optimal solution for fourth order systems requires 
continuous maximum control effort. For all but discrete sets of 
initial conditions, three switching points where control effort re- 
verses, are required. In this study, a quasi-optimal solution was 
obtained by using switching logic for only one switching point. The 
philosophy is to use maximum control effort during the initial 
portion of the solution and linear control for exponential settling 
during the terminal portion. Switching logic is used to terminate 
the maximum control effort at the proper time. This switch must 
take place sufficiently early to prevent excessive overshoot, yet 
late enough to result in a reasonably fast response. Control effort 
is allowed to decay during the terminal portion of the solution. This 
linear portion then replaces the last two switches of the minimum time 


solution and therefor prevents chatter mode, 
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In section 2.1 it was shown that the fourth order system could be 
represented in the uncoupled form as two second order systems. Ey 
doing this, the well developed second order switching lines of Titus 
and Demetry can be utilized to give an approximate type of switching 
logic for the fourth order system. 

Since two second order systems result from the uncoupling of the 
fourth order system, some sort of priority must be used to determine 
which uncoupled state pair should be applied to the second order switch- 
ing criteria. Various priority schemes were used including selection 
of the second order uncoupled state pair having higher velocity, higher 
energy, and variations of the latter. Also a switching Line which 
averaged the states of both the uncoupled state pairs was tried. 

The mathematical expression for minimum time second order switch- 
ing lines varies according to the type of plant eigenvalues involved 

(6 |. For this investigation, the second order switching line for 
null roots was used. This switching line is 


X,, | X 
control = -N sgn X Pl 
ON 


(ie 
oe) 
Cad 
See” 


fd 
1 
where N is the saturated control effort. 
A digital computer was used to investigate this type of control. 


See Appendix V for the digital computer program. In each of the 


following cases, the plant of Figure 1 with a = 0.5 and w = 1.0 wa 


> 


used, The terminal linear control was determined by the feedback 
coefficients described in equation (28) which give a dominant charac- 


teristic equation root of -0.2. 
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Figure 7 shows the system time response when only saturating 
linear control is used. Initial system positiom fer each curve is 
100 feet. Initial system velocity is varied from 30 feet per second 
to -30 feet per second. Control saturation occurs at plus or minus 
100. Time required for completion of 90 percent of the desired travel 
varies from approximately 7 seconds to 17 seconds for the extreme 
initial conditions. These travel times will be used to evaluate the 
following control systems which use full control during the initial 
portions of each solution. 

Figure 8 shows the system time response when the second order 
switching logic of equation (33) is applied to the average of the 
uncoupled positions and the average of the uncoupled velocities. 


The initial control equation therefore is changed to 


bee i Gre ce re 
control = -N sgn 1%; « Teel * Xu ial : 


(34) 
Initial conditions and control saturation remain the same. Time re- 
quired for completion of 90 percent of the desired travel varies from 
approximately 3 seconds to 10 seconds. This switching logic causes 

the system to respond approximately twice as fast as it did when only 
linear control was used, However, this control logic allows overshoot 
when the system has negative initial velocity. Therefore, it may be 
unacceptable for some applications. 

Figure 9 shows the system time response when the second order 
switching logic of equation (33) is applied to the uncoupled state pair 
having the higher velocity. Initial conditions and control saturation 
level remain the same. This shows that the use of the uncoupled 


velocity to determine the uncoupled state pair priority is not accept- 
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able. For initial system velocities of 10 feet per second or more, the 
response is identical to the linear response of Figure 7. For initial 

system velocities of zero feet per second or less, excessive overshoot 

occurs. 

Figure 10 shows the system time response when the second order 
switching logic of equation (33) is applied to the uncoupled state 
pair having the higher energy. To do this, each of the uncoupled 
velocity states was normalized by dividing by system natural frequency, 
w. Then the uncoupled state pair, whose states were farther from the 
phase plane origin, were assumed to have the higher energy. Initial 
conditions and control saturation level remain unchanged. A fairly 
consistent overshoot resulted from all the initial conditions con- 
sidered. Therefore, this control logic apparently has some potential. 
However, the logic must be altered to cause the switching point to 
occur sooner. 

Figure 11 is similar to that of Figure 10 except that the control 
logic was changed to cause the switching point to occur earlier, The 
initial control equation was changed to 

Yy |¥o| 


control = -N sgn Yay =a ee ; (35) 


A considerable improvement was realized by changing the control logic. 
The overshoot resulting from all initial conditions is now much 
smaller, Time required for completion of 9O percent of the desired 
travel varies from approximately 3 seconds to 5 seconds for the extreme 


initial conditions. 
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Figure 12 is similar to Figure 11 except the initial control logic 
was again changed to cause the switching point to occur earlier. The 


initial control equation was changed to 


Yo LY. | 
control = -N sgn yy +  OON : (36) 


Initial system positions are 100 feet, 70 feet, and 40 feet. Initial 
system velocities are 30, 0, and -30 feet per second, Moderate over- 
shoot occurs for small initial position combined with large negative 
velocity. Otherwise, the system is not extremely sensitive to small 
initial positions. This control system would therefore be acceptable 
for some applications where a small amount of overshoot could be 


accepted under extreme initial conditions. 
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4.0 Conclusions 

By selecting the proper closed loop system characteristic equa- 
tion roots, a fourth order system may be controlled to give a desired 
exponential response within initial condition limitations. The 
dominant characteristic equation root will completely define the 
system output trajectory after decay of the components associated 
with the non-dominant roots. Maximum system acceleration is a 
function of initial conditions and characteristic equation roots. 

A root locus study may be used to select characteristic equation 
roots which will minimize the possibility of oscillatory response 
components in the presence of state sensor errors. 

Within initial condition limitations, quasi-optimum time con- 
trol of a fourth order system is possible using second order switch- 
ing logic combined with terminal linear control. 

It is suggested that further studies be conducted to investigate 
the use of other types of switching logic for optimum control of 


fourth order system. 
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APPENDIX II 


Graphical Solution of System Acceleration 


Program SETROOT 


enter 











read 
initial 
conditions 


Calculate 
acceleration 





final 
time? 


© 


yes 


no F 
final 


curve? 


y 


yes 


call 
GRA PH 


(2) 
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APPENDIX III 


Gaussian Noise Generator 


Algorithm 
Given a uniform distribution p(x) = 1, -.55x-.5 


E(x) = uu = 0.0 


02 
Var(x) = , (x-u)°p(x dx = 1/12 


2d 
12 
Let R = aS x 
1 
E(R) = 0.0 
Le 
Var(R) = Var(cnx) = L Var( a x) 
> 1 
= Cc nh Var(x) 
2 
= C 
and 
p(R) = Normal(0,c“) by the Central Limit Theorem 


Computer Program Flow Chart 


sum 12 rand. 
numbers having 
Unit, dist. 

from -0.5 

LOLO25 


multiply 
by desired 
deviation 





ho 








Graph 1 Distribution of 1000 Random Numbers versus Theoretical 
Distribution with Standard Deviation of Unity 
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APPENDIX IV 


Linear Smoothing of Noisy Data 


Mathematical Development 
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minimize: y (y 
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‘e _ qc 
error 


define line y = Ax - B 


tf 


n 
2 2 
ober ag 2 ZY ~ Ax, - B) 


il 


(1) ae /dA = 2 20y cee ex) 


(2) fee /dB = x AM, = AK B)(-1) 


tl 


Ih 
>” n 
ye —ely, Z xX, = nx 
i ji= 
(2) B= y - Ax 
7. 73 
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xs -nxx 
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APPENDIX V 


Use of Second Order Switching Logic 
for Quasi-Optimum Control 


Program OPTFOUR 


enter 
read yes calculate 
initial linear 
conditions control 
has 
switch 
occurreg control yes 
saturated? 
no 
no 










calculate 
uncoupled 
states 
















calculate 
optimum 
control 






set 
: l 
inc flag 
time 









call 
RKUTTA 






Switch 


indicated? yes 






no 


no 


yes 


call 


‘ end 
GRA PH 
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